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We study the dynamics of a quantum superconducting circuit which consists of a Josephson 
charge qubit, coupled capacitively to a current biased Josephson junction. Under certain 
conditions, the eigenstates of the qubit and the junction become entangled. We obtain the 
time evolution of these states in the limit of weak coupling. Rabi oscillations occur, as a result 
of the spontaneous emission and re-absorption of a single oscillation quantum in the junction. 
We discuss a possible way to experimentally determine the quantum state of the junction and 
hence observe the Rabi oscillations. 

1 Introduction 

Advances in quantum information theory i are at the origin of the recent search for quantum 
mechanical two- level systems that can be used as quantum bits (qubits). In order to process 
quantum information in a useful way, certain operations on systems which consist of many qubits 
must be performed, such as the preparation and manipulation of, as well as a measurement on, 
entangled quantum states of coupled qubits. In this context, superconducting devices using small 
Josephson junctions are of particular interest. It has been experimentally demonstrated that a 
single Cooper pair box constitutes a two-level system which can be coherently controlled B'En-l. 
Moreover, the use of a Xlooper pair box as a qubit in the context of quantum computers was 
studied theoreticallyHfl'Q. But as of yet, the existence of entangled states, which are at the heart 
of quantum information processing, has not been demonstrated experimentally. 

In this article we study one of the simplest superconducting circuits in which entangled 
states can be realized!. It consists of a Cooper pair box coupled to a current-biased Josephson 
junction. Theoretically, this quantum circuit can be described by a two-level system coupled to 
a harmonic oscillator. After a description of the quantum circuit in the next Section, we discuss 
the time evolution of its eigenstates in Sec. ^. In particular, we demonstrate the existence of 
entanglement, very similar to the entanglement found for an atom in an electromagnetic cavityEI. 
In Sec. we discuss a possible way to measure the Rabi oscillations, associated with the dynamics 
of the entangled states. The last Section contains a discussion of the feasibility of our proposal. 

2 Superconducting circuit 

The circuit that we will study is depicted in Fig. |l[ It contains a quantum superconducting 
circuit connected to an electrodynamic environment. The quantum circuit itself consists of 
three elements, which we will discuss in some detail below. 

2. 1 Cooper pair box 

The Cooper pair box consists of a small superconducting island, coupled capacitively to a gate 
voltage Vg (gate capacitance Cg). The island is furthermore connected capacitively to a super- 
conducting electrode via a Josephson junction (capacitance Cj, Josephson energy Ej). Using 
the basis states |n), where n corresponds to the number of excess Cooper pairs on the island. 
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Figure 1: Superconducting quantum circuit connected to an electrodynamic environment. 



we can write the Hamiltonian for the box as 



i^box = Ec,j ^(2n - Ngf\n){n\ - ^ ^ (|n + l){n\ + \n- l){n\) . 



(1) 



Here, Ecj = e^/2Cj^eS is the total charging energy of the box (Cj^gg = Cj + [1/Cr + ^/{Cc + 
Cg)]~^) and Ng = —CgVg/e. The first term of this Hamiltonian corresponds to the electrostatic 
energy of the box. Generally, this energy is minimal if the number n of excess Cooper pairs 
on the island is an integer. However, if the dimensionless gate-voltagc Ng is an odd integer, 
Ng = 2m + 1, a degeneracy occurs, and the number n of Cooper pairs fluctuates between m and 
m + 1. The second term of the Hamiltonian, corresponding to the Josephson coupling between 
the island and the electrode, lifts this degeneracy. As a result, close to odd integer values of 
Ng, a gap opens up of order Ej and the relevant charge states are superpositions of the charge 
states \m) and \m + 1). In particular, if the gate-voltage is such that A^^ ~ 1, the states with 
n = and n = 1 are almost degenerate. At low energies, the Hamiltonian -f^box involves only 
these two states, and thus can be written as a matrix 



^^box ^ 



-Ej/2 Ec,j{2-Ngf 



The two eigenstates |— ) and |+) are superpositions of the charge states |0) and |1) 
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(2) 



(3) 
(4) 



where = 1 - /?2 = [i + 5Eg/^{5Egf + E]]/2 with 5Eg = -AEcjSNg and SNg = Ng - 1. 
The corresponding eigenenergies are 



E^ = Ec,j[l + {SNg)'] ^ -J{5EgY + £2 



(5) 



We see that the Cooper pair box behaves as a quantum mechanical two-level system, usually 
referred to as a Josephson charge qubit 



2. 2 Current biased Josephson junction 

A current biased Josephson junction can be characterized by two conjugate variables: the charge 
Q on the junction and the phase difference ^ across it. Ignoring effects related to the presence 



of the environment, we can write the Hamiltonian as 



Hr = ,^ + U{ct>), (6) 

where C^^efr is the effective capacitance of the junction in the circuit, C^^eS = Cr + [^/Cj + 
l/{Cc + Cg)]^^ , and U{(l)) = —Ej. cos(p— (<I>o/27r)/0 is the well-known tilted washboard potential 
($0 = h/2e is the superconducting flux quantum). It depends on the Josephson energy Er and 
the bias current /. For bias currents below the Josephson critical current Ic = 27r£'j./$0i the 
potential f7(</>) contains local periodic minima as a function of (j). These minima are quadratic 
with characteristic frequency ojr = {\/8E^^Ec^/h){l — {I /Ic)"^)^/^, where Ec^r = e^/2Cr,eff- 

For very low bias currents (I <^ Ic), the minima are separated by large potential barriers 
AU ~ 2Er — ^oI/2. As a result, in the quantum limit, various low-lying states can be found 
in each minimum if Er ^ Ec,r', the broadening of the energies of these states due to tunneling 
between neighboring minima can be ignored. Hence these states are well approximated by har- 
monic oscillator eigenstates, i.e., we can approximate Hr — hLOr{a^a + 1/2), where the operator 

(a) creates (destroys) an oscillation quantum in the junction. The eigenstates are denoted 
IXfc), corresponding to the presence of A; = 0, 1, 2, . . . oscillation quanta in the junction. Thus, 
at low bias current, the junction behaves as a quantum mechanical superconducting resonator. 
Since in this limit the phase is localized in a well-defined minimum of the potential U{(j)), the 
time-averaged voltage Vm across the junction remains zero. 

If the bias current is increased to values close to Ic such that Ic, the potential f/(</>) is well 
approximated by a cubic potential with a barrier height given by AU ~ (4\/2/3)£',.(l — ///c)^^^. 
Both AU and hojr decrease with increasing /; they vanish when / = Ic. The number of localized 
states in a given well decreases as AU/huJr ~ (1 — I/U)^^^. Moreover, the remaining levels 
broaden due to quantum tunneling out of the minima of the potential U {</)). A tunneling 
event changes the state of the junction0, thereby leading to a finite voltage Vm across it that 
can be detected. The broadening of the ground state energy is given by the tunneling rate 
Fq ~ tiLOrexp{—AU/hu)r). Since the excited states are located closer to the top of thebarrier, 
the tunneling rates increase with increasing level index k: Ffc/Fo ~ {432AU/hLOr)'''^^. The 
tunneling rates for the first two states \xo) and \xi) are therefore very different. We will see in 
Section ^ how this property can be used in order to distinguish these two quantum states. 



2.3 Coupling capacitance 

The capacitance Cc plays a crucial role in the circuit of Fig. |l| since it couples the charge qubit 
and the current-biased junction to each other. As a result, these two circuits are no longer 
independent and the system must be considered in its totality. Physically, the capacitance Cc 
couples the charge 2n — Ng on the qubit to the charge Q on the junction. Using the relation 

Q = ^fiuJrCr,cfi/2{a — a^)/i, the coupling Hamiltonian can be written as 

Hc = -iEconpi{2n-Ng){a^ -a). (7) 

Here we introduced the characteristic coupling energy Scoupi = \Jf>'^r/Ec,rEc,c/2, where Ec,c = 
e'^/2Cc,eS, with Cc,eS = Cc + Cg + (1/Cj + 1/Cr)~^. We will see in the next section how this 
coupling energy leads to entanglement between the states of the qubit and the junction. 



3 Operation at low bias current: entanglement 

In this section we discuss the behaviour of the circuit for values of the bias current away from 
the critical current, I <C Ic- In this limit, the system can be considered as a two-level system 



(the qubit) coupled to a harmonic oscillator {the current biased junction). The dynamics of 
such a system 1ms been discussed in the past 113. More recent applications are concerned with 
quantum optics □ or Josephson junctions 11^. Here we will show how entangled states can be 
obtained which involve both the qubit and the junction!. 

Throughout this section we will be interested in the situation Ng ~ 1, such that we have to 
consider the qubit states |— ) and |+) only. Furthermore, as far as the junction is concerned, we 
will limit the discussion to the ground state |xo) and the first excited state |xi)- Hence, in the 
absence of coupling, £^coupi = 0, the four lowest energy eigenstates are | — ,Xo)) \ — iXi)i l+iXo); 
and Let us now consider the special case fiuJr = Ej = E, i.e., the states \—,Xi) 

|+,Xo) are degenerate for Ng = 1. Their energy is E with respect to the energy of the ground 
state |— ,Xo)- In the limit of weak coupling, -Ecoupi E, this degeneracy is lifted: the states 
| — ,Xi) and |+,Xo) get mixed and hence give rise to the two entangled states and IV'2)) 

= [|-,Xi) + ^l+,Xo)]/^, 
IV'2) = [\-,Xi) -i\+,Xo)]/V2, 

corresponding to the energies E — Scoupi and E + -Ecouph respectively. 

Suppose that the system has been prepared in the state \tp{t = 0)) = |+,Xo) at time t = 0. 
This can be achieved by a suitable manipulation of the gate voltage Vg at times prior to t = oi. 
If we keep Vg fixed such that A'g = 1 at times t > 0, the time evolution of {ipit)) is given by 

|^(t)) = -L re-*(-e-Scoupi)t/?i|^^) _ e-i(^+£coupi)i/?i|^2)l . (8) 
\/2i L J 

We see that the state \ip{t)) oscillates coherently between |— ,xi) and |+,Xo)- In fact, these 
so-called quantum Rabi oscillations can be interpreted as the spontaneous emission and re-ab- 
sorption of one oscillation quantum by the junction. An interesting quantity is the probability 
Pxoit) to find the junction in the state |xo) (no oscillation quanta) after a certain time t. This 
probability shows Rabi oscillations as a function of t with frequency 2£'coupi/^) 

Pxoit) = \{+,xo\m)\' = \[l + cos{2E,oupit/h)]. (9) 

Since these Rabi oscillations are characteristic for the entanglement realized in the circuit, their 
measurement would provide direct evidence of the presence of the entangled states |'0i) and 
\ip2)- We will discuss the possibility to perform such a measurement in the next Section. 



4 Operation at bias current close to Ic'. quantum measurement 

In order to measure the probability Pxo{t) to find the junction in the state |xo) at time t, it 
is necessary to perform a series of measurements. Each of these measurements consists of the 
preparation of the state \ip{t = 0)) = |+, xo) at time t = and the subsequent coherent evolution 
until the time t, when the actual measurement of the state of the junction is performed. The 
statistics of the ensemble of measurements then yields Pxa{t). 

The proposed measurement procedure consists in transforming the two distinct quantum 
states Ixo) and |xi) into two different stable classical states which are easily measurable: the 
zero- voltage and the finite voltage state of the current-biased junction. This measurement is a 
one shot method which gives binary information on the initial quantum state of the junction. 

In order to measure the state of the junction at a certain time t, we propose to proceed as 
follows. Starting at time t, the bias current I is increased during a short ramping time St to a 
value Im close to Ic- The bias current is kept at Im for a time At and then switched back to a 
low value away from Ic- Finally, the DC-voltage across the junction is measured. 



The increase of the bias current to Im drastically increases the tunneling rates Tq and Ti 
of the two possible states of the junction. During the time At, the junction can transit into 
a finite voltage state as a result of a tunneling process. However, as we have seen, the escape 
time depends strongly on the initial quantum state. For instance, Fi ~ IOOOFq if AU/huJr > 3. 
Therefore, if the junction is in the state a DC voltage starts to develop after a time 1/Fi, 
whereas this will happen after a much longer time I/Fq ^ 1/Fi if the junction is in the state 
Ixo)- Hence if we switch the bias current / back to a value away from Ic after a time At such that 
1/Fi ^ At ^ I/Fq, the junction remains in the zero- voltage state if and only if the junction 
is in the state \xo)- It is during the time At that the two quantum states bifurcate into two 
different voltage states, hence we refer to this time as a measuring time. After At, the current 
is switched back to a lower value. Since the I — V characteristics of the junction are hysteretic, 
the zero voltage and finite voltage states remain dynamically stable. Thus we have a relatively 
long time to perform the actual voltage measurements and distinguish the two junction states. 

5 Discussion 

The functioning of the quantum circuit presented in this paper relies on various assumptions. 
Below we will discuss these assumptions for each part of the circuit in some detail. 

Josephson charge quhit. One of the main assumptions concerning the qubit is the absence 
of quasiparticles. If a quasiparticle tunnels onto the island, the state of the Cooper pair box 
changes: it no longer behaves as a quantum mechanical two-level system. Thus we need the 
superconducting gap to be sufficiently large, A > Ecj + Ej/2. We furthermore need the box to 
be quantum mechanically coherent during the experiment. This is one more reason to impose the 
absence of quasiparticles. But in addition it means that the coupling to the outside world must 
be sufficiently weak. Generally, the main source of decoherence is formed by time-dependent 
fluctuations of the effective gate charge Ng, either induced by fluctuations of the gate- voltage 
Vg, or by residual dynamics of charged impurities close to the qubit. However, if the qubit is 
operated close to the degeneracy point Ng = 1, the effect of this noise will be weak, since the 
energies E± depend only quadratically on the amplitude of the fluctuations. 

Current biased Josephson junction. When using the current biased Josephson junction 
as a quantum mechanical resonator, it is essential for the states |xo) and |xi) to be well- 
defined, i.e., the broadening Fo,i of the corresponding energy levels must be small compared 
to -Ecoupi- For this reason we need to work at low bias current, such that Fi ~ IOOOFq ~ 
1000^ ErEc,r exp ^ Er/Ec,r -^coupi- In addition, decoherence to the external circuit should 
not occur on a time scale of the order of h/Ecoupi- This means that the impedance should be 
sufficiently high. For a purely resistive environment for instance {R finite and L = in Fig. |l]), 
the typical relaxation time is RCr- We do require this time to exceed the Rabi period ?i/-Ecoupi- 

When using the junction as a detector it is important that it changes its state before re- 
laxation processes induced by the environment occur. Thus we need to increase the current I 
sufficiently close to Ic such that 1/Fi is much smaller than the relaxation time. 

Apart from increasing the tunneling rates F^, an increase of the bias current / has another 
important effect. Since the frequency ujr depends on /, the resonance condition huJr = Ej is 
no longer satisfied after the increase. As soon as the difference \huir — Ej\ exceeds ii^couph the 
qubit and the junction are no longer entangled and the states |— a-nd |-|-,Xo) become again 
eigenstates of the system. The probability to be in either of these states |xo) or ceases 
to evolve in time. This way, the probability P^.^ remains "frozen" into its value at the time t 
of the measurement (as long as relaxation phenomena can be neglected). In order for t to be 
well-defined, the ramping time 6t must be short compared to the period ?i/-Ecoupi of the Rabi 
oscillations. On the other hand, the ramping itself should not induce transitions in the junction, 
thus we need to impose the condition tOr^t 3> 1. 



In most experiments the environment is dominated by a resistor R. On the one hand, 
R should be large in order for the relaxation time to be long. On the other hand, R should 
be low enough to avoid heating at low temperature due to the bias current flowing through 
the resistor. In order to avoid these confliciting conditions, we propose to use a large on-chip 
inductor in series with the junction (see Fig. ||). To minimize the effect of the environment, 
the inductor has to be much larger than R/u>r and 1/(Ccj^). Then the relaxation time is 
given by {LuJr)'^C / R, which can be much larger than typical i?Cr-times t3; the characteristic 
frequency of the junction is not affected by the environment. Moreover, if the temperature T 
of the environment satisfies the condition R/L <^ 2kT/% <ti tOr, the escape rate is given by 
ln(ro(T)/ro) = {9/7rV2){h/e^){Cr/LIc){l-Im/Ic)^^'^kT, where To is the zero-temperature rate 
for quantum tunneling. The excess escape rate due to thermal fluctuations in the resistor is, 
surprisingly, independent of the resistance R of the environment. 

For the numerical estimates presented below we will consider parameters of typical super- 
conducting circuits li3'Q. For the Josephson charge qubit we have chosen Ej = 26.1/ieV and 
Ecj = 70//eV. The coupling capacitance is chosen to be of the order of Cj, Cc = 0.5fF, yielding 
Ec = 256neV. For a Josephson junction with a capacitance of 6.35pF and a critical current of 
10/iA, the characteristic frequency is about ISGHz at zero bias current. 

The resonant condition hcor = Ej for the charge qubit and the Josephson resonator is 
obtained at a bias current of 9.4^A. The period of the Rabi oscillations is given by T^abi ~8 ns. 
Taking into account the environment of Fig. [l| with R = 8$! and L=10nH, the relaxation time 
is about 136ns and thus longer than the period of the Rabi oscillations. 

During the measurement process, the bias current is chosen to be Im = 9.91//A, correspond- 
ing to a characteristic frequency of 4GIIz. The relaxation time is 51ns. For T = 20mK, we 
estimate the escape time to the voltage state to be 6//s from the state |xo) ^-nd 7ns from the 
state Ixi)- Taking At = 50ns, we expect the junction to bifurcate into the voltage state with a 
probability larger than 99.9% for the state and smaller than 0.1% for the state |xo)- 

In conclusion, we proposed an efficient procedure for single-shot quantum measurements of 
entangled states in a Cooper pair box coupled to a current-biased Josephson junction. 
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